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Abstract We present simulations of dwarf nova out- 
bursts taking into account realistic variations of the mass 
loss rate from the secondary. The mass transfer variation 
has been derived from 20 years of visual monitoring and 
from X-ray observations covering various accretion states 
of the discless cataclysmic variable AM Herculis. We find 
that the outburst behaviour of a fictitious dwarf nova with 
the same system parameters as AM Her is strongly influ- 
enced by these variations of the mass loss rate. Depending 
on the mass loss rate, the disc produces either long out- 
bursts, a cycle of one long outburst followed by two short 
outbursts, or only short outbursts. The course of the trans- 
fer rate dominates the shape of the outbursts because the 
mass accreted during an outburst cycle roughly equals the 
mass transferred from the secondary over the outburst in- 
terval. Only for less than 10% of the simulated time, when 
the mass transfer rate is nearly constant, the disc is in a 
quasi-stationary state during which it periodically repeats 
the same cycle of outbursts. Consequently, assuming that 
the secondary stars in non-magnetic CV's do not differ 
from those in magnetic ones, our simulation indicates that 
probably all dwarf novae are rarely in a stationary state 
and are constantly adjusting to the prevailing value of the 
mass transfer rate from the secondary. 

Key words: accretion, accretion discs - binaries: close - 
stars; individual: AM Her - novae, cataclysmic variables. 



1. Introduction 

Cataclysmic variable systems (CV's), in which a white 
dwarf accretes material from a Roche lobe filling sec- 
ondary star (see Warner 1995 for an enyclopaedic re- 
view) can be broadly divided into three subclasses: in 
non-magnetic systems, the white dwarf accretes via an 
accretion disc; in the magnetic polars, the infalling mat- 
ter couples directly onto the strong magnetic field of the 
white dwarf before it can build a disc and is funnelled 
to accretion region(s) near the magnetic pole(s) of the 
white dwarf; finally, in CV's containing a weekly mag- 
netic white dwarf (intermediate polars), a partial disc may 



exist, with the mass flowing from the inner edge of the 
disrupted disc through magnetically funnelled accretion 
curtains onto the white dwarf. 

Dwarf novae are non-magnetic cataclysmic variable 
stars which show characteristic eruptions with photomet- 
ric amplitudes in the range of 2-8 magnitudes Q These 
eruptions typically last a few days to several weeks and 
recur quasi-periodically on timescales of weeks up to many 
years. They are generally thought to result from thermal 
instabilities associated with hydrogen ionization in the ac- 
cretion disc (see Cannizzo 1993a for a review and Ludwig 
et al. 1994 for a detailed parameter study). 

An important boundary condition of the disc-limit- 
cycle model is the mass transfer rate from the secondary 
to the accretion disc. In most studies of dwarf nova out- 
bursts this mass transfer rate is kept constant (Cannizzo 
1993b; Ludwig & Meyer 1997; Hameury et al. 1998, 1999). 
Duschl and Livio (1989) were the first to examine com- 
bined mass transfer and disc outbursts, though within the 
framework of the mass transfer instability model where 
individual and short-lived mass transfer events are capa- 
ble of producing single outbursts. Smak (1991, 1999) dis- 
cussed mass transfer variations in the context of the su- 
peroutburst phenomenon in dwarf novae of the SU UMa 
type and dwarf nova outbursts with enhanced mass trans- 
fer during outburst. 

King & Cannizzo (1998) and Leach et al. (1999) tested 
how the accretion disc in a dwarf nova system behaves if 
the mass transfer from the secondary varies abruptly be- 
tween different levels. They found that these mass transfer 
variations produce only subtle effects on normal dwarf no- 
vae, including variations in the outburst shape, and that 
the dwarf novae keep on having outbursts even if the trans- 
fer rate is reduced close to zero. 

In spite of these efforts we are left with the question of 
what the real variations of the mass tranfer rates from the 
secondaries in dwarf nova systems are and how they can 
influence the outburst behaviours of the accretion discs. 
Fortunately, nature provides an answer to the first ques- 



^ Also intermediate polars are known to show outbursts in 
their partial discs; e.g. EX Hya, GK Per 
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tion in the form of discless cataclysmic variables, the po- 
lars or AM Her systems. In these systems, the mass trans- 
fer rate can be estimated directly from the observations 
because there is no accretion disc acting as a mass buffer. 
Thus, we can use the long-term light curve of an AM Her 
system as a measure for the mass transfer variations in a 
fictitious but realistic dwarf nova system with the same 
system parameters. In the present paper, we present the 
results of such a numerical experiment. 

In the following section, we briefly review the equa- 
tions for the viscous and thermal evolution of the accretion 
disc, discuss the numerical methods used in our code, and 
compare the results of two standard models with other 
fine mesh computations. Thereafter we derive the mass 
transfer rate in AM Her as a function of time, Mtr(i)- Fi- 
nally, we apply this mass transfer rate to a fictitious dwarf 
nova and discuss the effects that can be observed in the 
outburst behaviour. 



Wc solve Eqs. (1) and (2) using a combined Finite- 
Element / Finite-Difference algorithm (FE for the spatial 
part and FD for the time-evolution). Apart from our own 
work (Schreiber & Hcssman 1998) the method of Finite El- 
ements has not been used in this context. As this method 
proved to be extremely robust, it warrants a somewhat 
more detailed description. 

The idea of FE is to divide the region of interest 
(the disc radii between i?in and .Rout) into n — 1 ele- 
ments and to expand the function u{x) which is supposed 
to solve the differential equation with suitable functions 
"(^) = SiLi o,iiPi{x) for every element. In order to get a 
continuous solution over all elements, the functions {ip) of 
every element have to be transformed to the so-called lo- 
cal basis functions Ni and the coefficients to the so-called 
nodes [ci] before collected together (Gruber & Rappaz 
1985, Schwarz 1991). To solve the differential equation, 
the function 



2. Finite Element Methods in the context of disc 
evolution 

The classical equation describing the viscous evolution of 
the surface density S in a geometrically thin, axisymmet- 
ric accretion disc is obtained by combining the vertically 
averaged Navier-Stokes and mass-conservation equations: 



+ 



2nR 



(1) 



where jy = {2/3)aRgT/ fifl is the kinematic viscosity, a 
denotes the viscosity parameter, Rg the gas constant and 
H the mean molecular weight. The first term on the right 
hand side is the classical disc diffusion term and the sec- 
ond term describes the external mass-deposition. For the 
purposes of this study, we follow Cannizzo (1993b) and 
simply add the mass lost from the secondary Afir(t) in 
a Gaussian distribution at a fixed radius near the outer 
edge. 

Similarly, one can derive an energy equation (in the 

central temperature T: 



dT 2(H-C + J) RJT 1 dlRvYi) dT 



dR' 



(2) 



dt CpE jjLc-p R dR 

where Cp denotes the specific heat, 

—^vd\og{vY.r^) 
'"^ " ~R d\ogR 

is the local radial flow velocity, H = represents 
viscous heating, C = crT^g is the radiative cooling and 

3 S a dT 
= r-'^RdR^'^dR^ 

the radial energy flux carried by viscous processes (see 
Smak 1984; Mineshige & Osaki 1983; Mineshige 1986; 
Ichikawa & Osaki 1992 and Cannizzo 1993b, for discus- 
sions). 



(3) 



fc=i 



has to meet the requirement formulated by Garlerkin: 

the integral of the residuum (which one gets by inserting 
Eq. (3) into the differential equation) weighted with the 
functions Nj{x){j = l,...,n) has to vanish. This require- 
ment, the interchange of integration and summation and 
partial integration lead to matrix-equations of the form 



Bc + Ac = D, 



(4) 



with A = {ttij), B = bij and D = di {i = 1, .., n;j = 1, .., n 
for n nodes). 

To solve the differential Eqs. (1) and (2), we have to fill 
A, B, Din the sense mentioned above and calculate c from 
Eq. (4). After transforming to the variables X = 2 J? 2 and 
S = XT, we derive for the surface density from Eq. (1): 



n 

fe=l 



+ 



bij = j NiNj 



X2 dX " dX 

12 dN.dN 

idX 



d. = Y.l '-^^N.dX. 



fc=i 
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Similarly it is easy to obtain for the central temperature 
Eq. (2): 

^ / f (i)dN, dN, 

fc=i 

+ / Pk^Nk^NjdR + j pfNiNkNjdR) 
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bij = I NiNjdR 
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Figure 1. Comparison with Cannizzo (1993b, his Fig.6d). 
From top to bottom we used 40, 60, 80, 100 and 200 nodes. 
Convergence is obtained for 100 nodes. 



k=l 



pfNuNjdR. 



The coefBcients p^^^ are given by 



(4) 



IJCpR dR 



off 



CpE 



The index k refers to the value of the quantities at node 
number k, which is equivalent to the radius Rk- Notice 
that we only used linear basis functions in this paper. 

As mentioned above, Eq. (4) has to be solved to get 
c{t + At) from c{t). Using simple finite differences leads to 

c{t + M) = c{t)]^M{B-'^ Ac{t) 



(5) 



In order to test our code, we carried out two sets of 
calculations using the binary parameters and cooling func- 
tions from Cannizzo (1993b) 



0.02, 

4.0 X 10^°cm, 



Ml — IMq, Qfh = 0.1, a, 
i?i„ = 5.0 X lO^cm, i?out 
A/tr = 1.5 X lO^Mg/yr 

and Ludwig & Meyer (1997) 

Ml = O.63M0, an = 0.2, = 0.04, 

i?i„ = 8.4 X lO^cni, i?out = 1-7 x 10^°cm, 

Aft,, = 5 X lO^^g/s. 

The resulting light curves are shown in Figs.j^ and ^ Our 
code reproduces the sequence of only relatively long out- 
bursts found by Cannizzo (1993b) for the parameter of 



Figure 2. Comparison with Ludwig & Meyer (1997). 
From top to bottom we again used 40, 60, 80 ,100 and 
200 nodes. Our code acurately reproduces the sequence of 
one long outburst followed by two short outbursts. 

SS Cygni as well as the sequence of one long outburst fol- 
lowed by two short outbursts found by Ludwig & Meyer 
(1997) to describe VW Hydri. The short outbursts arise 
when there is not enough mass stored in the disc and there- 
fore the heating wave gets reflected before it has reached 
the outer edge of the disc. 

We find that at least 100 nodes are necessary for long- 
term convergence. The outburst and quiescence duration 
decreases with an increasing number of nodes because - 
with finer zoning - the length of time spent on the viscous 
plateau becomes shorter (Cannizzo 1993b). This effect is 
smaller in Fig. |2| because the disc is smaller in this system. 

We conclude that our FE-code produces results which 
are in excellent agreement with those of other fine-mesh 
computations. 

3. The mass loss rate of the secondary star in 
AM Herculis 

The strong magnetic field of the white dwarf primary in 
polars prevents the formation of an accretion disc. With- 
out an accretion disc acting as a buffer for the transferred 
mass, the mass loss rate from the secondary equals the 
mass accretion rate on the white dwarf at every moment, 
Afacc = AJftr (the free-fall time is ^ Ih). As an obser- 
vational consequence, any variation in rate at which the 
secondary star loses mass through the Li point will re- 
sult in a quasi-immediate change of the observed accre- 
tion luminosity. The brightest polar, AM Her, has been 
intensely monitored at optical wavelengths by observers 
of the AAVSO for more than 20 years and shows an ir- 
regular long-term variability, switching back and forth be- 
tween high- and low states of accretion on timescales of 
days to months. 

The problem of deriving the mass loss history of the 
secondary star is then equivalent to that of determining 
the accretion luminosity Lace ^ a function of time. As the 
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Figures. The mass transfer rate in AM Her as a function of time. 



bulk of the accretion luminosity is emitted in the X-ray 
regime, a bolometric correction relating the densly moni- 
tored optical magnitude to the total luminosity has to be 
derived. This approach has been followed in detail by Hess- 
man et al. (A&A, submitted) using X-ray observations ob- 
tained at multiple epochs. We summarize only briefly the 
results here. The accretion luminosity is computed from 
the observed accretion-induced flux F„rr as 



iacc(<) =47rd2Facc(i) 



(6) 



with d ~ 90pc (Gansicke et al. 1995). We further use 
Face ~ ^sx + 3 X Fhx with Fsx and f hx the observed soft 
and hard X-ray fluxes respectively. The factor three ac- 
counts for the additional cyclotron radiation emitted from 
the accretion column and for the thermal reprocession of 
bremsstrahlung and cyclotron radiation intercepted by the 
white dwarf and emitted in the ultraviolet (Gansicke et al. 
1995). The mass loss ( = transfer) rate is then 



J-'acc (t)i?wd 



(7) 



where G is the gravitational constant and i?wd and M^d 
are the white dwarf radius and mass, respectively. As the 
actual properties of the white dwarf in AM Her are still 
the subject of controversial discussions (Gansicke et al. 
1998, Cropper et al. 1999), we use the parameters of an 
average white dwarf, M^d = 0.6 and i?wd = 8.4 x 
10* cm. Mtr(t) is shown in Fig.^. The average value of the 
mass transfer rate in AM Her is Mav = 7.88 x 10-^^gs"^ = 
1.24 X lO~^°M0yr~^. The derived mass transfer rates of 
AM Her are in general agreement with results published in 
the literature. For example, Beuermann & Burwitz (1995) 
found transfer rates between 0.8 and 2.0 x 10~^° AfQyr"^ 
and Greeley et al. (1999) estimated a mass transfer rate 
of 2 x lO^^gs"^ for the high state of AM Her from far 
ultraviolet spectra. 



4. Results 

4-.1. The fictitious dwarf nova 

We devised a fictitious dwarf nova with a non-magnetic 
primary of mass M^d = 0.6 Mq and an orbital period of 
P — 3.08 hr, i.e. a non-magnetic twin of AM Her. For these 
binary parameters, we obtain Rm — Rwd = 8.4 x 10* cm 
for the inner disc radius and Rout = 2.2 x 10^" cm for the 
outer edge of the disc. We then used our FE code to follow 
the structure of the accretion disc in our fictitious dwarf 
nova for 7000 d, applying the variable mass transfer rate 
Mtr{t) derived above and standard viscosity parameters 
ah = 0.2 and ac = 0.04. 

We show 500 day-long samples of our calculations in 
Figs. 4-8. In each figure, the top panel show the mass 
transfer rate as a function of time (solid line) and the av- 
erage mass transfer rate log(Mav[gs~^]) = 15.90 (dotted 
line). The panel below displays the disc mass Mdisc nor- 
malized with the averaged disc mass Mdisc — 1-64 x lO^^g. 
The two lower panels display the light curves calculated 
with the varying mass transfer rate and the light curves 
calculated with the constant average mass transfer rate 
respectively. For the constant average mass transfer rate 
the disc goes through a 60 day-long cycle including one 
long outburst followed by two short outbursts. The long 
outbursts are those in which the entire disc is transformed 
into the hot state while the short outbursts arise when the 
outward moving heating wave is reflected as a cooling wave 
before it has reached the outer edge of the disc. 

Our numerical experiment clearly demonstrates that 
the outburst light curve of the flctitious system is strongly 
affected by the variations of the mass transfer rate. Even 
in the case of relatively small fluctuations effects on the 
outburst behaviour (16.0 < log(Mtr[g s~^]) < 16.4) are 
clearly present in the light curve. In Fig.|^, the mass trans- 
fer rate is always high during the 500 days but the disc 
switches between an accretion state with only long out- 
bursts and states where one or two short outbursts follow 
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Figure 4. From top to bottom as a function of time: 
the mass transfer rate (top, solid hne), the averaged mass 
transfer rate (top, dotted hne), the normahzed disc mass 
and the hght curves produced by the fictitious dwarf nova 
with the variable transfer rate and the averaged transfer 
rate adopted from AM Her. 



a long outburst. When the tranfer rate decreases some- 
what (~ day 200), the disc does not save enough mass to 
create consecutive long outbursts until the mass transfer 
rate increases again (~ day 350). 

In addition to this effect, our experiment shows that a 
sharp decrease in the mass transfer rate instantaneously 
changes the outburst behaviour of the accretion disc. In 
Fig. H the disc first behaves as in Fig. ^ but when the trans- 
fer rate drops sharply (day 4340), the long outbursts im- 
mediately vanish and the duration of the quiescent phase 
increases somewhat. 

Another remarkable point is that even during rel- 
atively long periods of very low transfer rates (Fig.^, 
days 4700-5000) the disc does not stop its outburst ac- 
tivity but produces only short outbursts with slowly de- 
creasing amplitudes and increasing quiescence intervals. 
This confirms the findings of King & Cannizzo (1998). 

Fig. 7 shows 500 days of our simulation in which the 
adopted mass transfer from the secondary varies strongly 
on a timescale of roughly 20 days. This is the most fre- 
quent case during our calculation but there are rare peri- 
ods of nearly constant mass transfer. Fig. 8 shows the light 
curve of the fictitious system during 500 days in which the 
mass transfer rate nearly exactly equals its averaged value 
(top). As a result the light curves computed with the real 
mass transfer rate and the averaged value look equal. 

In summary, one can say that the variations of the 
mass transfer rate leads the disc to switch between three 
states in which only long outbursts occur (log(Aftr) > 
16.3), one long outburst is followed by one or two short 
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Figure 5. The same as Fig. ^ but another snapshot of the 
simulation. The sharp decline of the mass transfer rate is 
immediately reproduced by the accretion disc. 
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Figure 6. The same as Fig.^ but another snapshot of 
the simulation. Even the long period of low transfer rates 
does not stop the outburst activity of the disc. 



outbursts (16.3 > log(Mtr) > 15.7), and only short out- 
bursts occur (log(Mtr) < 15.7). 

In order to understand the described behaviour of the 
disc we take into account the viscous timescale ty ^ B? jv 
which gives an estimate of the timescale for a disc annulus 
to move a radial distance i?. For the quiescent state tv,c 
and the outburst state tv,h and with R = i?out (where the 
mass transferred from the secondary is added to the disc) 
we obtain for the viscous timescale 

<vx 2000 d, (8) 

tv,h 15 d. (9) 
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The mass added to the disc during quiecence is stored in 
the disc because it moves inward on the long timescale tv,c 
whereas during an outburst even mass from outer regions 
can reach the white dwarf within a viscous timescale. This 
makes it possible that the mass accreted onto the white 
dwarf during a long outburst can be up to roughly one 
third of the disc mass day 4220). Therefore the disc can 
relax to equillibrium with the mass transfer rate in only 
one outburst in the case of high transfer rates (Fig. 4) . 

The prompt response of the disc to the sharp decline 
in the mass transfer rate (Fig. 5) can be understood in the 
same way: due to the short viscous time (iv,h) the disc ac- 
cretes a substantial fraction of the disc mass {r-^ l/4Mdisc) 
during the last long outburst which immediately prevents 
long outbursts when the mass transfer rate becomes low 
(~ day 4340). 

Finally, the long period of low transfer rates (Fig. 6) do 
not prevent outbursts because the mass accreted during 
the short outbursts is only a few percent of the disc mass. 
Therefore the disc relaxes to the mass transfer rate on a 
longer timescale. 

4-. 2. The mass transfer and mass accretion rates 

An important point in understanding the physics of ac- 
creting binaries is to know how far the outburst behaviour 
and hence the resulting light curves depend on real varia- 
tions of the mass tranfer rate. 

To answer this question, we compare the averaged 
mass accretion rate onto the white dwarf with the mass 
transfer rate. Figure 9 shows that the time in which the 
disc relaxes to an equilibrium with the mass transfer rate 
depends on the occurence of long outbursts: when the 
accretion rate is averaged over 20 days (dotted line in 
Fig. 9) the mass transfer and the accretion rate correspond 
roughly only during periods where only long outbursts 
occur (days 40 -150, see also Fig. 4) but for the periods 
where the disc goes through a cycle of short and long out- 
bursts the accretion rate has to be averaged over 60 days 
to match the mass transfer rate (days 150-500). If the 
mass transfer rate drops steeply and stays in a low-state 
(days 4700 -5000 in Fig. 10, see also Fig. 6), the accretion 
rate needs more than 60 days to follow this behaviour, 
because the disc produces only short outbursts in which 
only a small percentage of the disc mass is involved. 

In order to make this plausible we give a relaxation 
timescale as the ratio of the viscous timescale in out- 
burst with the relative mass fraction accreted during an 
outburst: 



tr 



t. 



(10) 



For high mass transfer rates this timescale is around 
Todays and so the correspondence of the averaged ac- 
cretion rate and the mass transfer rate in Fig. 9. is not 
surprising. In the case of low transfer rates (Fig. 10, day 
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Figure 7. The same as Fig.^ but a snapshot of the sim- 
ulation where the transfer rate strongly varies. 




Figure 8. The same as Fig.^ but a snapshot of the sim- 
ulation where the transfer rate nearly stays constant. 



4700-5000) this timescale is longer (~ 300 days) because 
only roughly 5 percent of the disc mass are accreted during 
an outburst. 



4-.3. Dependence on the primary mass 

In the numerical experiment above, we have assumed an 
average white dwarf mass for the primary in AM Her. The 
literature holds a large spectrum of white dwarf mass es- 
timates for AM Her, M^d = 0.39 Mq (Young & Schneider 
1981), Mwd = 0.69 M0 (Wu et al.l995), M^d = 0.75 Mq 
(Mukai & Charles 1987), M^d = 0.91 Mq (Mouchet 1993) 
and Mwd = 1.22 Mq (Cropper et al. 1998). Based on 
the observed ultraviolet spectrum of AM Her and on its 
well-established distance, Gansicke et al. (1998) estimated 
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Figure 9. The mass transfer rate (solid line) and the 
accretion rate onto the white dwarf as function of time. 
The accretion rate is averaged over 60 days (dashed line) 
and averaged over 20 days (dotted line). 
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Figure 10. The mass transfer rate (solid line) and the 
accretion rate onto the white dwarf as function of time. 
The accretion rate is averaged over 60 days (dashed line) . 



i?wd ~ 1-1 X 10^ cm, and, using the Hamada-Salpeter 
(1961) mass-radius relation for carbon cores, 0.35 Mq < 
-^wd < 0.53 Mq. As the Hamada-Salpeter mass-radius 
relation is valid for cold white dwarfs, the finite temper- 
ature « 20 000 K of the white dwarf in AM Her would al- 
low also somewhat higher masses, M^d ~ 0.65Mq, which 
is very close to the average mass of field white dwarfs, 
0.6Mq, that we used. 

Even though we exclude a massive white dwarf based 
on the observational evidences, we repeated our simulation 
with Mwd = 1.0A/q in order to test the sensitivity of our 
results. 

In a first step, we recompute Mtr(i) from Eq. (|^) with 
-Mwd = 1.0 Mq and a corresponding i?wd = 5.4 x 10^ cm. 
The resulting mean accretion rate is 3.0 x lO^^gs^^, a fac- 
tor 2.6 lower than before. Then, we simulated once more 



0.1, Qfc = 0.02 and i?out = 2.8 x 10^ 

In Fig.|Tl| we show 500 days of our calculation with 
M^vd — I.QMq. The disc produces only short outbursts 
and the outburst cycle of four outbursts with decreasing 
amplitude is hardly changed even by drastic variations of 
the mass transfer rate (day 4700). 

The different responses that our fictitious dwarf novae 
with 0.6Mq and I.OMq white dwarfs show to the vari- 
able mass transfer rate are easy to understand: both the 
increased primary mass and the decreased radius of the 
white dwarf reduce - as mentioned above - the derived 
average mass transfer rate. 

In addition, the outer disc radius of the fictitious sys- 
tem with Mwd = I.OMq increases. The disc becomes more 
massive, i.e. Mdisc — 4.59 x lO^'^g, because of its increased 
size. 

Due to the reduced mass transfer rate and the in- 
creased disc size, the heating waves are able to reach the 
outer edge of the disc only during the first outburst of the 
outburst cycle and only in cases where the mass transfer 
rate is high (log Mtr > 16.1). Even in this rare situation 
the disc stays only a few days in the hot state and ac- 
cretes only a small fraction of the disc mass ('^1/8 Mdisc) ■ 
For lower transfer rates (the extremly more frequent case 
shown in Fig. 11) the heating front gets reflected before 
it has reached the outer edge of the disc. Hence, only a 
small percentage of the disc mass is involved in any out- 
burst. Therefore, and due to the longer viscous timescale, 
the relaxation timescale given in Eq. (10) is always larger 
than a few years. 

Summing up, the adopted primary mass and the aver- 
age accretion rate play an important role on the influence 
that the variable mass transfer rate has on the outburst 
behaviour. 



5. A note on the disc limit-cycle model 

A number of new features have been added to the well 
known limit-cycle model (Meyer & Meyer-Hofmeister 
1981; Smak 1982; Cannizzo, Gosh & Wheeler 1982). 
Hameury et al. (1998, 1999, 2000) have shown how the 
light curves change if the disc size is allowed to vary and 
which effects irradiation has on the outburst behaviour. 
Additionally these authors discussed other sources of un- 
certainties such as the tidal torque and the evaporation 
of the inner parts of the disc (Meyer & Meyer-Hofmeister 
1994). They concluded that all these effects must be in- 
cluded in order to obtain meaningful physical informa- 
tion on e.g. the viscosity from the comparison of pre- 
dicted and observed light curves. Another point of im- 
portance is the interaction of the accretion stream leav- 
ing the secondary star and the accretion disc. Schreiber & 
Hessman (1998) tested the influence of stream overflow on 
the disc evolution in dwarf novae. They found that signif- 
icant stream overflow can lead to reversion of the inward- 
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travelling cooling front and create an outward-travelling 
heating front. This behaviour would produce small dips in 
the light curve during the declining phase. 

The model we used here contains neither irradiation 
nor allows it the disc radi to vary. The stream mass is 
simply deposited in a small Gaussian distribution near 
the outer edge of the disc. Therefore we do not attempt 
at present a comparison with observed light curves. Nev- 
ertheless, all calculated light curves with constant mass 
transfer from the secondary relax in a quasi-stationary 
outburst cycle (of one or more outbursts) which repeat 
periodically. Our model shows in a qualitative way that 
real mass transfer variations may have a dominant influ- 
ence on the outburst behaviour at least in systems with 
relatively small discs and strong mass transfer. 

6. Discussion and conclusions 

There is no reason to assume that the mass transfer varia- 
tions of the secondary observed in AM Her are not present 
in non-magnetic systems. Our numerical experiment in- 
cluding realistic variations of the mass transfer rate in a 
dwarf nova system is, therefore, a significant step towards 
a better understanding of dwarf nova light curves, and, 
thereby, of the underlying disc limit cycle. 

The light curve produced by our fictitious system 
switches between three states depending on the actual 
mass transfer rate. High transfer rates lead to only long 
outbursts where the entire disc is transformed into the 
hot state. If the transfer rate is near the average value 
the disc goes through a cycle of three outbursts, one long 
outburst followed by two short ones. Even long periods of 
low transfer rates do not force the disc to stop its outburst 
activity: long outbursts are suppressed and the duration 
of quiescence increases but the disc always produces short 
outbursts. From this follows that the low-states of VY 
Sculptoris stars (a subgroup of novalike variables) could 
not be caused by low transfer rates alone (see also Leach 
et al. 1999). 

We find that in our fictitious system the mass accreted 
during an outburst cycle is dominated by the course of 
the mass transfer rate if the mass transfer rate varies 
significantly. The disc always relaxes to equilibrium with 
the mass input from the secondary. Thus, our experiment 
strongly supports King & Cannizzo's (1998) claim that 
dwarf nova accretion discs are probably never in a sta- 
tionary state but are constantly adjusting to the prevailing 
value of Mil-. Only during periods where the mass transfer 
rate is nearly exactly constant the disc periodically repeat 
the quasy stationary outburst cycle. Such periods are rare 
but occur in AM Her. 

The strong influence of the mass transfer rate on the 
outburst behaviour of the fictitious system clearly indi- 
cates that probably most (if not all) the deviations from 
periodic outburst cycles seen in the light curves of dwarf 
novae are caused by variations of the mass transfer rate. 
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Figure 11. The same as in Fig.g but assuming a more 
massive primary M„d — 1.0 Mq 
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